Jet correlation measurement in heavy-ion collisions: 
from RHIC to LHC 



O . Constantinos A. Loizides 



o 

(N 



X 



Massachusetts Institute of Technology, 77 Mass Ave, Cambridge, 02139, MA, USA 
E-mail: loizides@mit.edu 



Oh' 

(D , 

^0 , Abstract. We attempt to deduce simple options of 'jet quenching' phenomena in heavy-ion collisions at 

. -y/sNN = 5.5 TeV at the LHC from the present knowledge of leading-hadron suppression at RHIC energies. 

' In light of the nuclear modification factor for leading particles we introduce the nuclear modification factor 
for jets, -R^A' f^"^ longitudinal momenta of particles along the jet axis, -R^a- 
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^ ! 1. Introduction 

^ I At RHIC energies, hard probes, mainly light quarks and gluons, are experimentally accessible in heavy- 

lO ' ion collisions with sufficiently high rates. Their production is quantified using inclusive single-particle 

O I spectra at high momentum (pT > 2 — 3 GeV). The suppression of these high-momentum (leading) 

^ ' particles in central nucleus-nucleus relative to peripheral or pp collisions is one of the major discoveries 
at RHIC. In Au-Au collisions at various centre-of-mass energies {^/s^ = 62 A, 130 and 200 GeV) the 

{iLiI two experiments with high transverse-momentum capabilities, PHENIX and STAR, but also PHOBOS 

^ ■ and BRAHMS, have measured: 

• the suppression of single-particle yields at high ( ~ 2 — 4 GeV) [1-5]; 

• the disappearance, in central collisions, of jet-like correlations in the azimuthally-opposite side 
(away-side) of a high-pT leading particle [6, 7] and, quite recently, the reappearance of the particles 

I on the away-side manifested in low-momentum hadrons [8, 9]; 

The absence of these effects in d-Au collisions at the same centre-of-mass energy {^/s^ = 200 GeV) [1, 
10-12] confirms that in (central) nucleus-nucleus collisions final-state (as opposed to initial-state) effects 
modify the measured particle spectra. 

The experimental observations have been explained in terms of various quenching models, where the 
energetic partons produced in the initial hard scattering 'lose' energy as a consequence of the interaction 
with the dense, partonic matter. Most models implement parton energy loss according to medium- 
induced gluon radiation of a hard parton traversing dense partonic matter of finite size [13-21]. Also 
hadronic interactions [22, 23] have been investigated and partially found to contribute to the observed 
depletion of the hadron spectra. 

At the LHC at 30 times higher centre-of-mass energy, hard probes will be abundantly produced, even 
at energies of more than one order of magnitude higher than at RHIC. These energetic partons might 
be identified by their fragmentation into hadronic jets of high energy. In contrast with RHIC, the initial 
energy of E't > 50 GeV is high enough to allow the full reconstruction of the hadronic jet, even in 
the underlying heavy-ion environment [24]. Therefore, measurements of changes in the properties of 
identified jets in nucleus-nucleus with respect to pp collisions will become possible. 



In the present work we attempt to deduce perspectives for 'jet quenching' measurements at LHC 
given present data and models of hadron suppression at RHIC. In section 2 we introduce PQM, a model 
that describes high-px suppression at RHIC, and state its prediction of the nuclear modification factor 
for LHC. In section 3 we argue that the surface-emission effect limits the sensitivity of leading particles 
to the density of the medium. In section 4 we present a PYTHIA simulation of jet quenching at LHC 
conditions and introduce the nuclear modification factor for identified jets, i^^A' ^^'^ ^'^^ longitudinal 
momenta of particles along the reconstructed jet axis, -R^A- 



2. Leading-particle spectroscopy at RHIC and LHC conditions 

The effect of the medium on the production of a hard probe is typically quantified via the nuclear 
modification factor, 

1 d'^Nl^'^/dpT di] 
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which measures the deviation of the nucleus-nucleus from the superposition of independent 
nucleon-nucleon collisions. In absence of strong nuclear- initial-state effects it should be unity, if binary 
collision scaling holds. 

Figure 1 (left) shows Raa{pt:) for central Au-Au collisions at y^s^ = 200 GeV as published by 
PHENIX [3, 4] and STAR [2] together with a pQCD calculation of paiton energy loss obtained with our 
Monte Carlo program PQM (Paiton Quenching Model) [20]. 
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Figure 1. Left: Raa{pt) for central Au-Au collisions at ^/i^ = 200 GeV. PHENIX [3, 4] and 
STAR [2] data ai^e reported with combined statistical and pT-dependent systematic errors (bars on the 
data points) and pT -independent systematic en^ors (bars at Raa = 1)- The theoretical calculation is 
obtained with PQM resulting in an average transport coefficient of 14 GeV^/fm. The band denotes 
the systematical uncertainty of the model. Right: The PQM prediction for Raa{pt) in central Pb-Pb 
collisions at y^s^ = 5.5 TeV. The LHC prediction presented in Ref. [17] (Vitev-Gyulassy) and PQM 
for central Au-Au collisions at RHIC, ^/s^ = 200 GeV, are shown for comparison. 



PQM combines a recent calculation of the BDMPS-Z-SW paiton energy loss [25] with a realistic 
description of the collision geometry at mid-rapidity. Our approach allows to calculate the transverse 
momentum and centrality dependence of single- and di-hadron correlation suppression, as well as the 
'energy-loss induced' azimuthal anisotropy of particle production in non-central collisions. The model 
has one single parameter, k, that sets the scale of the density of the medium. Once the parameter is 
fixed, e.g. on the basis of the nuclear modification data at ^/s^ = 200 GeV, we scale it to different 
energies and collision systems assuming its proportionality to the expected volume-density of gluons. 
The application of our model is restricted to the high-px region, above 4 — 5 GeV at RHIC and above 
10 GeV at LHC energies, since we do not include initial-state effects. 



The leading-particle suppression in nucleus-nucleus collisions is obtained by evaluating 
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in a Monte Carlo approach. Equation (2) describes the production of high-px hadrons within the 
perturbative QCD collinear factorization framework, at mid-rapidity, r] = 0, including medium-induced 
parton energy loss in the BDMPS-Z-SW formalism, /^(fe) is the parton distribution function for a parton 
of type a{b) carrying the momentum fraction rr^^;,), is the partonic hard-scattering cross sections 

and Di^ij{zj) is the fragmentation function, i.e. the probability distribution for the parton j to fragment 
into a hadron h with transverse momentum pT = Zjptj. The modification with respect to the pp 
(vacuum) case is given by the energy-loss probability, P{AEj; Rj,ujcj), for the parton j. Its two input 
parameters, the kinematical constraint and the characteristic scale of the radiated gluons, depend on the 
in-medium path length L of the parton and on the BDMPS transport coefficient of the medium, q. The 
latter is defined as the average medium-induced transverse momentum squared transferred to the parton 
per unit path length. (See Ref. [20] on how we calculate the two parameters R and lOc for a given 
parton and realistic geometry.) In the original calculation [25], the energy-loss probability (quenching 
weight) is independent of the energy of the original parton allowing a parton with finite energy to radiate 
more than its energy. To account for finite paiton energies we have introduced two different ways of 
constraining the weights, the non-reweighted and reweighted case. In the first case we constrain the loss 
to the energy of the parton, whenever the radiated energy is determined to be larger than that. In the 
second case we require by truncation of the distribution that the energy loss is less than the energy of the 
paiton. The resulting energy loss is larger in the non-reweighted case, because partons 'thermalize' with 
a probability de P{e). It is argued [21, 25] that the difference in the values of the observables for the 
two approaches illustrates the theoretical uncertainties of the model. 

Coming back to fig. 1 (left) we note that for the chosen value of k the calculation reproduces Raa{pt) 
for central Au-Au collisions at y/s^ = 200 GeV. In Ref. [20] we consistently compare model 
predictions to various high-pT observables at RHIC energies. We also show that we need to scale k 
by a factor of seven to obtain the LHC prediction of i?AA at ^/s^ = 5.5 TeV shown in fig. 1 (right). 
Our prediction for the LHC is in agreement, both in the numerical value and in the pT-dependence, with 
that obtained in Ref. [21], while it is quite different from the one calculated in Ref. [17]. 



3. Properties of leading-hadron suppression 

At LHC energies, the expected Raa for central Pb-Pb collisions at LHC as a function of px is rather 
flat, i.e. almost px -independent, similar to what is observed at RHIC at px > 4 GeV. To study this px- 
dependence we show in fig. 2 Raa{pt) in central Pb-Pb at LHC for (q) = 12 GeV^/fm (corresponding 
to the value found at RHIC), (q) = 24 GeV^/fm and (q) = 98 GeV^/fm (the value obtained by the 
scaling from RHIC to LHC), as well as the result of a calculation with fixed q = 10 GeV^/fm and fixed 
length of 4.4 fm. Clearly, the fixed-length calculation shows a stronger px-dependence than the PQM 
calculation. 

As we explain in the following, the flatness of Raa is a consequence of 

/ \ "(pt) 

• the steeply falling production cross-section, oc ( ^^eihts^ j > where n(px) is rising from about 7 
to 12 (RHIC) and from 6 to 7 (LHC) in the relevant px regime; 

• the emission from the surface, which for large medium densities dominates [26]. 
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Figure 2. i?AA as a function of px for 0-10% most central collisions at LHC energy obtained by PQM. 
The calculations in the parton-by-parton approach (solid lines) are compared to a calculation for fixed 
transport coefficient and length (dashed). All graphs are in the non-reweighted case. 



Raa, eq. (1), at mid-rapidity, can be approximated by 



^aa(pt) = I dAEP{AE, PT + AE) 



dNPP{pT + AE) I diVPP(pT) 



dp- 



dpi 



(3) 



where, dA^P^/dpT is the spectrum of hadrons (or partons) in the case of no medium (i.e. pp neglecting 
initial state effects). The suppression computed with eq. (3) is found to give a rather good approximation 
to the one computed with PQM or with the full formula, eq. (1). In the case the production 
spectrum is (approximately) exponential the px-dependence cancels in the ratio and we find Rkk to 
be (approximately) independent of pT- At RHIC this is the case at about px > 30 GeV. Below that 
value and at the measured values of 5-12 GeV, as well as at the LHC (n(pT) < 7), the spectrum is given 
by a power-law. It was expected [27] that TiAA oc (1 + c / \/ ?^(pt) Pt)""^*''^^ i-C- reaching unity for 
large px. 

However, this neglects the fact that for dense media surface emission or, more generally, the 
probability to have no energy-loss, P{AE = 0, E), plays a significant role, an effect which is even 
more pronounced at low px (compai^ed to uJc)- To simplify our argumentation we allow either no loss 
{AE = 0) or complete loss (AE = E) in the non-reweighted case, P{AE, E) = po 5{AE) + (1 - 
Po) 6{AE — E). ^ Inserting the constrained weight into eq. (3) we obtain 



7?aa(pt) =P* + (1-P* 



diVPP(2pT) / diVPP(pT) 



dpi 



dpi 



(4) 



It is obvious that eq. (4) is just a crude approximation, but it demonstrates that the value of i?AA is 
dominated by the fraction of hadrons (or partons), which escape without losing much of their energy. 
For the simple power-law production spectrum the contribution from higher px is suppressed by about 
(1 + AS/pt)"^^^-*. Taking into account only fixed values of q and L the probability p* is given by the 
discrete weight, po, at 72 = 0.5 q L^, amounting to about 0.05(0.002) for quarks(gluons). However, for a 



^ For a dense medium the constrained weights at low parton energy indeed do have a sharp peak at zero and at maximum 
possible energy loss, whereas the values in between are negligible. 




Figure 3. Distributions of parton production points in the transverse plane (upper row) and in-medium 
path length (lower row) for partons that escape the medium and produce hadrons with pT > 50 GeV in 
central Pb-Pb collisions at 5.5 TeV for {q) = 12 GeVVfm (left) and {q) = 98 GeVVfm (right). The 
quantity (L) denotes the average of the path-length distribution. The calculation is for the non-reweighted 
case. 
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Figure 4. Raa as a function of the average transport coefficient, (q), for 10 and 100 GeV hadrons in 
0-10% central Pb-Pb collisions at ^/s^ = 5.5 TeV. The calculation is for the non-reweighted case. 



proper calculation one must take into account the right production ratio of quarks-to-gluons. For realistic 
path-length distributions p* is dominated by partons, which are emitted close to the surface and, thus, 
enhanced relative to po obtained at fixed scale. It turns out that p* evaluated with PQM at LHC central 



conditions, averaged over path-lengths and parton types, is independent of in the range shown above, 
and takes values of about 0.14, 0.1 and 0.05 for the scales used in fig. 2. 

Surface emission will be present even for very large-momentum hadrons. Figure 3 (top) visualizes 
the distribution of production points (xq, yo) in the transverse plane for partons, which escape from the 
dense overlap region and yield hadrons with pl^^'^^°^^ > 50 GeV. The chosen values of the medium 
density correspond to the value found at RHIC and the expectation for LHC. The corresponding in- 
medium path-length distributions, fig. 3 (bottom), reveal that the average 'thickness' of the emission 
surface is restricted to less than 2 fm: even hadrons with > 50 GeV ai^e emitted dominantly from the 
surface. 

The dominance of the surface effect limits the sensitivity of leading particles to the density of the 
medium, mainly for experimentally accessible low-px range at RHIC. This is demonstrated in fig. 4 
where we show the dependence of Raa on the average transport coefficient evaluated with PQM for 
0-10% most central collisions at ^Snn = 5.5 TeV first shown in Ref. [21]. For 10 GeV hadrons the 
nuclear modification factor is sensitive to average medium densities up to about 15 GeV^/fm, and for 
100 GeV hadrons the sensitive regime might widen to average values of about 30. 

4. Jet spectroscopy and modification of jet properties 

At LHC, leading-hadron spectroscopy will be naturally extended by jet spectroscopy. Simulations show 
that jets with transverse energies of more than 50 GeV are identifiable and reconstructible on an event- 
by-event basis, even in most central collisions; however with severe limitations in the resolution of the 
jet energy. The interaction of the hard partons with the medium created in these collisions is expected to 
manifest in the modification of jet properties deviating from known fragmentation processes in vacuum. 
Calculations predict that the additional gluons radiated by the original parton remain (fragment) inside 
the jet cone, although redistributed in transverse phase space [25, 28]. The corresponding jet-production 
cross section is expected to follow binary scaling. However, the jet shape is claimed to broaden and the jet 
multiplicity to soften and increase [29]. Ideally, if one reconstructs the hadronic energy for outstanding 
high-energy jets the jet energy may be associated with the original energy of the parton. 

In combination with results at lower parton energies from coiTclation methods jet measurements at 
higher energy might complete the picture of medium-induced parton-energy-loss phenomena. In the 
following it is our aim to introduce a simulation and analysis of jet quenching effects. Details can 
be found in Ref. [24]. Unquenched jets representing the pp reference are generated using the PYTHIA 
generator [30]. Quenched signal jets ai^e simulated as an afterburner to PYTHIA developed by A. Morsch. 
It introduces parton energy loss via final-state gluon radiation, in a rather ad-hoc way: Before the partons 
originating from the hard 2-to-2 process (and the gluons originating from ISR) aie subject to fragment, 
they are replaced according 

partonj(^) partoni(^ - A^) + n{AE) g\uon{AE/n{AE)) 

conserving energy and momentum. For each parton, AE is calculated by PQM in the non-reweighted 
case and depends on the medium density, paiton type and parton production point and emission angle 
in the collision overlap region. The number of radiated gluons, 1 < n{AE) < 6, is determined by the 
condition that each gluon must have less energy than the quenched parton from which it was radiated 
away. First results using this afterburner at LHC energies have been discussed in Ref. [31]. 

For the following sample analysis we have prepai^ed events containing jets simulated with the modified 
PYTHIA version for (q) = 1.2, 12 and 24 GeV^/fm, which are embedded into 0-10% central 
HIJING [32] events. The choice of first value implies only a very little modification of the embedded 
quenched jets compared to embedding of jets with standard PYTHIA. The second corresponds to the 
case found to describe the Raa at RHIC, whereas the third is a conservative choice about a factor of four 
smaller than the expected value extrapolated from RHIC to LHC. 

The jet finding is based on final particles (no detectors effects) at mid-rapidity (— 1 < < 1). We 
use the standard ILCA cone finder [33]. However, to cope with the soft heavy -ion background we need 



restrict the cone radius to R = 0.3. For charged particles we require pt > 2 GeV. The same settings 
are used for the reference jet measurement in pp. In central Pb-Pb collisions the mean reconstructed jet 
energy is about 80% for 50 GeV jets, and more than 90% for 100 GeV jets; the resolution is about 20% 
and 10%, respectively. 



4.1. Nuclear modification factor for jets 

It has been shown [29] that the medium-induced broadening of the jet reduces the energy inside R = 0.3 
by 15% and by ~ 7% for jets with E't = 50 and E't = 100 GeV, respectively, and already at 
R = 0.7 the effect reduces to about 2%. However, one must be cautious about these findings, since this 
prediction has been calculated assuming a rather low value of the gluon density at LHC conditions. 
We have verified (without embedding in HIJING) that for i? = 0.7 more than 90% of the energy 
remains in the jet cone. Up to 45% of the rise at jet energies lower than 50 GeV can be attributed 
to the contribution of the underlying heavy-ion background to the jet signal. Jet suppression (at fixed 
value of R) may be observed in the suppression of inclusive single-jet spectra in central nucleus-nucleus 
relative to peripheral or pp collisions. Analogous to eq. (1), we therefore define the nuclear modification 
factor for jets, R^^j^{ET, rj; h), according to 



R^i{ET,v; b) 



1 
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(5) 



Equation (5) implicitly depends on the (fixed) value of R used to find the jets. Figure 5 shows R-j^j^{ET) 
for different quenching scenarios in 0-10% central Pb-Pb. 
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Figure 5. The nuclear modification factor for jets, R''j^j^{Et), as a function of the reconstructed jet 
energy, E^'^, for different quenching scenarios in 0-10% central Pb-Pb. The jets in pp and Pb-Pb ai^e 
identified with the cone finder using R = 0.3 and pT > 2 GeV. (Full calorimetry is assumed.) 



4.2. Longitudinal momentum modification factor 

The signature of medium-induced gluon radiation should be visible in the modification of the jet 
fragmentation function as measured through the longitudinal and transverse momentum distributions of 
associated hadrons within the jet. The momenta pai^allel to the jet axis, pL = Phadron cos(^jct, ^hadron), 
ai^e expected to be reduced (jet quenching), while the momenta in the transverse direction, jx = 



Phadron sin(^jct, 6'hadron). to be increased (transverse heating). In our simulation for technical reasons 
we can not give an explicit transverse-momentum kick to the radiated gluons. Therefore, we concentrate 
on the changes of the longitudinal fragmentation. 

We define the longitudinal momentum modification factor according to 



(6) 



where for proper normalization A^^A ^pp^ account the total number of jets that were found in 
A-A and pp, respectively. Equation (6) implicitly depends on the transverse jet-energy range of jets 
taking into account in the distributions. Figure 6 shows i?AA(PL) for jets with reconstructed energies of 
Ej?"^ > 80 GeV in 0-10% central Pb-Pb compared for the different quenching scenarios. 
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Figure 6. Longitudinal momentum modification factor, -R^A' for jets with Ej^'^ > 80 GeV and different 
quenching scenarios in 0-10% central Pb-Pb. The jets in pp and Pb-Pb are identified with the cone 
finder using R = 0.3 and > 2 GeV. (Full calorimetry is assumed.) 

The expected behavior is clearly visible: higher medium density leads to stronger suppression for 
large longtitudinal momenta and enhancing of smaller momenta. At low pL the effect becomes becomes 
most apparent. At these low momenta there is an moderate additional contribution from particles which 
belong to the underlying heavy-ion background, which amounts to about 20% for pl < 5 GeV. It might 
be interesting to study the (average) slope of -RaaIpl) as a function of the reconstructed jet energy. 

5. Conclusions 

We have presented the PQM model [20] which is able to describe high-px suppression effects at RHIC. 
Its application to LHC conditions leads to the expectation that Raa is essentially constant with p^, 
also for very high momenta up to 100 GeV. The nuclear modification factor for leading hadrons is 
largely dominated by surface effects. This limits its sensitivity to the density of the created medium. The 
measurement of reconstructed jets above the underlying heavy-ion background may allow to probe the 
medium to deeper extents. In a simulation of jet quenching at LHC we have studied two observables: 
the nuclear modification factor for jets, -RaA' ^^'^ the nuclear modification factor for the longitudinal 
momenta of particles along the jet axis, i?AA- "^^^ definitions implicitly depend on the parameters used 
to find the jets. For jets with E't ^ 50 GeV both ai^e suitable measures to quantify deviations from 
the pp case: R-'j^j^{ET) decreases with increasing E^ and i?AA is enhanced at low and suppressed 



at high PL- However, for jets with Et < 50 GeV suffer from the contamination of the jet cone 
by uncorrelated particles from the underlying background which at these energies severely influences 
measured jet properties. 
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